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We develop a framework for the multiple scattering of a polarized wave. We consider 
particles with spin propagating in a medium filled with scatterers. We write the amplitudes 
of each spin eigenstate in a local, mobile frame. One of the axes is in the direction of propa- 
gation of the particle. We use this representation to define a directional Green's operator of 
the homogeneous medium and also to write the spin-dependent scattering amplitudes. We 
show that this representation reveals a Berry phase. We establish a generalized Green-Dyson 
equation for the multiple scattering problem in this framework. We show that the general- 
ized Green-Dyson equation can be solved by linear algebra if one uses a representation of 
the rotations based on Wigner D-matrices. The properties of light scattering are retrieved if 
we use spin 1 particles. Our theory allows to take into account several kinds of anisotropics 
like circular or linear dichroism and birefringence, Faraday effects and Mie scattering within 
the same formalism. Several anisotropics can be present at the same time. 



I. INTRODUCTION 

In optics, several phenomena involve the polarization of light. Such phenomena, like birefrin- 
gence, dichroism and Faraday effect, are well understood and described in textbooks yj], usually 
in the case where the medium is homogeneous and linear. When filled with particles, the medium 
is no longer homogeneous and in addition to the transport of light, one has to study its scattering 
by the particles. A century ago, Mie discovered that the scattering of light by small dielectric 
spheres is linear and also affects the polarization . The theory of light scattering has been 
extended for several different particles and studied in the context of multiple scattering [|4|,|5|,|6|]. 
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While the polarization is well studied for a single scattering, the problem becomes more difficult 
in multiple scattering and a scalar approximation is often used. Such an approximation makes it 
simpler to investigate the coherent effects using sophisticated diagram techniques 0,180. As a con- 
sequence, polarization effects are neglected whereas they may, in certains situations, be of interest. 
In magneto-chiral birefringence, for instance, a residual degree of circular polarization persists is 
chiral media [SJ] . Recent experiments of coherent backscattering of light by cold atoms exhibit 
polarization-dependent enhancement factors [10]; this result is interpreted thanks to the quantum, 
nonlinear, nature of the cold atoms [JJJ]. Under certain conditions, it was even forecast that an 
anti-localization effect could be observed Jl^ . Quantum effects of light were also introduced, in 
second quantization, for the study of dense cold atomic gases jl3\\ . 

Remarkably, most of the scalar approximation techniques are applicable not only to light but 
also to microwaves, ultrasound and non interacting electrons, as emphasized in Ref. OD. The exis- 
tence of a polarization, alike, is not restricted to light; other waves display polarization effects: The 
other electromagnetic waves, elastic waves, electrons. More generally every particle with a spin is, 
from a quantum point of view, a polarized wave. The classical light polarization is a consequence 
of the vectorial nature of the electromagnetic field whereas in the quantum-mechanical picture, 
photons have two states called helicities. The origin of the helicity of a particle follows from the 
spin. Relativity imposes that the spin state of massless particles does not exist in the direction 
of propagation, photons therefore posses only two helicities. The freedom of the relative phase 
between the two helicity states is at the origin of the phenomenon of polarization as we experience 
it in optics. All quantum particles with a non-zero spin also have helicities and hence may display 
polarization effects. Therefore, the problem of polarization transport in multiple scattering does 
not only arise in optics: Polarization transport in multiple scattering is a more general problem and 
optics can be seen as the most investigated particular case. 

According to the quantum duality of particles and waves, we may describe the propagation 
of particles as wave amplitudes from a statistical, or probabilistic, point of view. Reciprocally, 
the propagation of a wave can be described by the statistics of particle paths, through techniques 
like Feynman's path integration. Following this picture, many theoretical calculations related 
to the multiple scattering of polarized waves are performed by Monte-Carlo techniques. Such 
techniques have been developped for light, even in the presence of quantum scatterers and 
for elastic waves Ill5ll . Monte-Carlo techniques are a pragmatic solution to the problem of taking 
into account each trajectory as a whole. Is it possible to consider each trajectory entirely and 
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perform statistics in a more formal way ? A pioneering work on trajectory statistics was made by 
Sato |Ql with the introduction of a directional Green's function for the elastic energy (equivalent 
to the optical intensity) in a three dimensional medium with nonisotropic scatterers. However, as 
his study concerned energy, polarization was not included in the model. 

The trajectories of energy, or intensities, have been studied in a multiple scattering context 
using Stokes parameters. These parameters are preferentially used by experimentalists because 
they are, as intensities, directly measurable quantities. They depend on the choice of the reference 
frame used for their observation. Their ability to take coherence effects into account is limited to 
intensity correlations. In order to dispose of the reference frame indermination, Kuscer and Ribaric 
introduced harmonical functions to describe the Stokes parameters Jl7|l . Based on Wigner's work 
concerning group theory [18], these functions have also been used under the name "spin- weighted 
spherical harmonics", though with a different phase factor II 1911 . They are also implicitely used in 
light scattering studies under the name "vector spherical harmonics" [3, 5I. The seminal work of 
Kuscer and Ribaric motivated a large number of studies on the phase matrix [Q, 



the radiative transfer equation Q23 



21, 



22D and on 



240. 



The role of trajectories appears in a surprising effect involving polarized light: A four-foiled 
pattern is observed in backscattering experiments of polarized light from a medium containing 



anisotropic scatterers II25LI261 127I1. In these experiments, the source is not a planar, infinite wave, 



but a localized beam. The effect does not exist if the source is extended on a scale of the order of 
the four-foiled patterns. The same medium does not exhibit such patterns if the observations are 
made in transmission rather than in reflection, so that the properties of polarization depend not only 
on the detailed trajectory, but also on the direction of observation, or more precisely on the frame 
of observation. It has been suggested that the appearance of these patterns is due to the existence 
of a Berry phase for the phot on |28I1 . In quantum mechanics, Berry phases have been studied 
already in many situations [29, 30]. In optics, the literature is abundant in experimental [31] and 
theoretical considerations [ 32 , 33!] . If light is not guided in an optical fiber, but rather multiply 



scattered, the Berry phase still exists Q3_4{| . The probability distribution of the phase depends on 
the statistics of the paths followed by the photons. 

The purpose of the present article is to introduce a general framework for the multiple scattering 
of any polarized wave, as simply and as pedagogically as possible. We will base our construction 
on the concept of path statistics, in the spirit of the Monte-Carlo calculations, or Feynman's path 
integral, to account for local interactions and path-dependent effects. We express the formalism 
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for particles with an arbitrary spin S, although we consider only linear media and scatterers and 
thus leave the other quantum properties to further investigations. Our objective is not to provide 
a complete theory for multiple scattering of polarized waves but rather to set up a framework in 
which such a theory could be developed. On that account we assume the simplest possible point 
of view and make several assumptions, that are briefly discussed. 

We have organized our work as follows: The representation of the field and the geometry of 
rotations are discussed in Section HB A directional Green's operator is discussed in Section (TTTl 
it describes the transport in a homogeneous medium. Using the same formalism, we introduce 
the scattering operator between incoming and scattered state in Section [IV] These first three sec- 
tions constitute the elements of the theory. We demonstrate that our theory takes into account the 
Berry phase without requiring any extra tool in Section |V} We derive a generalized Green-Dyson 
equation for multiple scattering of polarized waves in Section [VI] and we show that a solution of 
this equation is obtained by means of the rotational harmonics, introduced in Section IVIIl We 
discuss the particular case of rotational invariant scatterers in a rotational invariant medium in sec- 
tion IVIIII and conclude. Fundamental formulae concerning the rotational harmonics are provided 
in appendix |Al and some of the several notations used for them are listed in appendix El 

H. GEOMETRY AND POLARIZATION 

We begin the presentation of our work with the introduction of the geometry elements used 
throughout this article. Mainly, the novelty of our work relies on the representation of the state of 
waves in local frames. Polarization is introduced as a transformation property of the wave state 
under local rotations. We express all space coordinates in a fixed reference frame XYZ. Any local 
frame X'Y'Z' is therefore obtained after a rotation of the reference frame. The coordinates of the 
local frame unit vectors, written in the reference frame, form a 3 x 3 matrix, which the unique 
rotation matrix mapping XYZ into X'Y'Z 1 and preserving orientation. The group of rotations is 
SO (3). Considering the rotation associated to a frame is mathematically equivalent to considering 
the frame itself. We henceforth only use the description in terms of matrices for mathematical 
convenience and use exclusively the ZYZ Euler representation of rotation matrices, such that any 
rotation R of SO (3) is decomposed into 



R = Z(0)Y(0)Z(V), 



(1) 
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where Z(0) is the rotation around the Z-axis of angle <fi. Similarly Y is a rotation around the 
F-axis. The decomposition (QQ) is unique if < 9 < n and < 0, ip < 2n. For 9 = or n, the 
decomposition is not unique. 

In our discussion, we use the concept of trajectory, or path of a particle. Studying the propaga- 
tion of a wave using the particle picture would not be fully general if we were concerned by the 
behavior of single particles. But as our goal is to describe all possible trajectories as a whole, it is 
known that this way of addressing the problem of multiple scattering is correct. As a consequence, 
we will simultaneously use the wave picture, as equivalent to the statistics of particle trajectories, 
and the particle picture of the propagating wave. The spin of the particle emerges as an essential 
ingredient of the theory, after geometry considerations. 

We study a medium filled with N scatterers and write for the position of the scatterer i. 
The trajectory of the particle is a succession of displacements between points i-j and changes of 
direction at the points r,. We call such a change of direction a scattering event. We consider the 
momentum p n after n scattering events and decompose it into 

p n = p n P„ Z (2) 

where p n = |p n | is the modulus of p„, P„ is a rotation matrix (expressed in the reference frame) 
and z is the unit vector of the reference frame along the third coordinate. We choose the local 
frame for the segment of the path after the n th scattering event and before the next one. The third 
coordinate of the local frame is therefore always pointing in the direction of motion. 

We remark immediately that the decomposition © of p„ is not unique: The Euler angle ij) of P n 
can be taken arbitrarily, it is the spin gauge freedom. Here stands a crucial point in our model: The 
introduction of the local frame is a natural way of introducing the spin, and hence the polarization, 
into a multiple scattering formulation. If the wave has a spin S, there are 2S + 1 possible values for 
the spin component along the direction of propagation, that we label with "s" (— S < s < S). At 
a given point r, the field is a superposition of partial fields, which we consider to be plane waves, 
with different momenta. For simplicity, we may consider all directions of propagation for a fixed 
absolute value of the momentum. The amplitude of probability for observing the spin state s in the 
frame P at point r and time t is written ^(r, t, P, s). In this article, we regularly use the bracket 
notation for functions: The functions / of a variable x are written as (x \ f), instead of f(x). There 
is no consequence of this notation concerning the physics itself, but we find it more convenient for 
our purpose, in particular for the introduction of the rotational harmonics in section [VO We then 
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write the amplitude as a bracket product 



v|/(r, t, P, s 



) = (r, P, s 



(3) 



Polarization is related to the frame of observation and if this frame is changed, polarization is 
modified according to certain rules. Light polarization, for instance, is turned by an angle — ip 
if the frame of observation is turned by an angle ip around the direction of propagation. It is 
therefore natural to change the basis of representation and use the eigenbasis of rotation along the 
direction of propagation. Light has two circular polarization states, corresponding to the helicities 
of the photon. A rotation of angle tp creates a phase shift for each circular eigenstate. If s is 
an eigenvalue of the spin operator along the direction of propagation, rotations along this axis 
commute with the spin operator. Unitarity implies that 



on the on-shell states. After these remarks, we shall name the third Euler angle, usually noted ip, 
the spin angle. 

We have presented a new description of the field at a given point that depends on the local 
frame and thereby taking into account the spin. We still have to check the completeness of the ket 
representation |P, s). As we have remarked, the direction of propagation is Pz, so that |P, s) and 
|P', s') are orthogonal if Pz 7^ P'z. From relation (0Q) we conclude that non orthogonal |P, s) and 
|P', s') have the same Euler angles except the spin angle, which is arbitrary. Let R and R' be two 
reference frames, we introduce the function <5j_(R _1 R') to denote 



where 9 and (f> are the first two Euler angles of R and 6' and </>' the first two Euler angles of R'. 
S± is a kind of Dirac function for rotation matrices, which select frames having the same third 
unit vector. Using this notation, the condition Pz = P'z is imposed by defining (P', s'\P, s) 
proportional to <5 ± (P _1 P'). We remark that 5 ± (P /_1 P) = <5 ± (P _1 P'). Finally, the orthogonality of 
the spin eigenstates yields (s\s'} = 5 SS >. From the expression © we get the product 



(PZ(V0, s 




(4) 



<J ± (R _1 R') = 4tt5(cos0' - cos 9)5(<j)' - 0) 



(5) 



<P', s' 



P, s) = e 




(6) 



where ip and ip' are the spin angles of P and P' respectively. 
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The superposition of the partial fields for all directions of propagations is expressed by 




dP v|/( r , t, P, s)e 




/ siaBde / d0 / V(r, t, Z(0)Y(0)Z(^), s)e 



SO (3) 



— / sin0d0 / d0M/( 
4Wo Jo 



r, t, Z(0)Y(0), s) ee v|/ s ( r , t). 



(7) 



The presence of the term e 1S1 ^ is essential. Its role is to add the amplitudes in such a way that the 
amplitudes in two distinct frames with the same direction of propagation are in phase and do not 
cancel out. Forgetting this term would make all ty s vanish except ty , because the exponential term 
from © would be integrated to give 0. Multiple scattering theories for scalar waves actually only 
consider the term ty . The completeness of the Hilbert basis |P, s) is expressed by the closure 
formula (1 is the identity operator) 



The integral over SO (3) has been defined in Eq. ©. The description we have introduced in this 
section will be used in the formulation of the multiple scattering for polarized waves. In the next 
section we introduce the Green's operator for the polarized states |P, s) and its space dependence 
in the absence of scatterers. 



Before studying multiple scattering, it is necessary to investigate the transport in a homoge- 
neous medium. By the word transport, we mean the response at a position r' and time t' in the 
frame P' (see previous section for the introduction of the frames) to a source at position r and time 
t in the frame P. It is advantageous to use the Green's functions because the two situations of mas- 
sive and massless particles can both be handled with the same formalism. In the case of a massless 
particle, like the photon, transport follows from Maxwell's equations and is often formulated by 
the Helmholtz equation. Massive particles transport is of a different nature and the dynamics of 
their wavefunction responds to the Schrodinger equation. 

We denote by Q the Green' function of the operator + ^ or of the operator — |^ Jp- — E 
to describe the transport of the electromagnetic field or of the wavefunction respectively. The 




P, s) (P, s\ = 1. 



(8) 



III. TRANSPORT IN A HOMOGENEOUS MEDIUM 
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equivalence between these two operators is obtained through the relation 

- = -V2mE, (9) 
c n 

so that we shall use in our formulae only the "massless" notation uj/c. The results we present can 
be extended to massive particles of arbitrary spin by using equation © and changing the sign of Q. 
If the medium is invariant under translation and in time Q(r, t; r', t') depends on r'-r and t' — t. 
In this section, we construct a Green's operator G which depends on the direction of motion at r 
and r', and we call it the directional free Green's operator. 

We introduce the free Green's operator as the operator transforming the wave function along the 
propagation of the wave if no scattering events occur. The transition amplitudes characterizing the 
response may depend on the spin. If we note G the free Green's operator, transport is described 
by the matrix elements (r ; , P', s'\ G (t', t) |r, P, s). We have introduced the position ket |r) and 
denoted |r, P, s) = |r)®|P, s). Naturally, it is not physical in quantum mechanics to consider the 
position and the direction of motion of a particle simultaneously. The directional Green's operator 
that we need for our theory can be seen as an intermediate element of computation: The physical 
Green's function is the superposition of the directional Green's operators for all initial and all final 
directions of motion. In Feynman's picture of path integrals, it corresponds simply to decompose 
the path integral formulation of Q into path integrals over trajectories with a constraint on the 
direction of motion at the initial and final point. We have illustrated this feature in figured] 

In a medium with translational invariance, the momentum is conserved so that Go has a factor 
(P|P'). Between two scattering events, the wave travels in space from r to r' along the direction 
p = Pz. The direction of the momentum, Pz is the same as the direction r' — r because of the 
relation r' — r = c(t' — t)Pz, so that we have to impose this constraint to the free Green's operator. 
To express the directional constraint on r' — r, we write a similar equation as © : 

r' - r = rDz (10) 

so that the dependence of the Green's operator G on D is simply reduced to 5j^(DP _1 ). Finally, 
we have the expression for the dynamics Green's operator 

(r', P', s'\ G dyn (t', t) |r, P, s) = £(r, t; r', f) ^(DP" 1 ) (P', s'\P, s) . (11) 

The factor S± (DP -1 ) in expression (fTT|) is a very strong restriction imposed by rotational invari- 
ance. It states that a rotation of the medium must be accompanied by the same rotation of the 
frames P and P' (through the term (P', s'|P, s)) to leave the Green's operator unchanged. 
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FIG. 1: Schematic view of the directional Green's operator. Without scatterers (upper figures) the only 
situation for which the Green's operator does not vanish is when the three vectors r' — r, Pz and P'z are 
colinear and have the same direction (l.h.s. of the figure). In this situation, only one path is possible: The 
straight line from r to r' sketched as a dashed line on the r.h.s. of the figure. 

In a complex medium (lower figures) the directional Green's operator (r', P'| Go |r, P) describes the trans- 
port from r to r' with the constraint that the initial and final directions of motion are Pz and P'z respectively. 
There are several different paths, some examples of which are shown as dashed lines on the r.h.s in the fig- 
ure. Note that the vectors r' — r, Pz and P'z need not to be colinear anymore. 



In the presence of scatterers, it has been observed in backscattering configurations the transport 
of polarization depends on the relative directions of the incident polarization beam and the vec- 
tor r' — r Q26Q - Moreover, in some cases, like in the presence of linear birefringence or dichroism, 
the direction of propagation has a strong influence on the transport. Such complex geometrical 
dependences will be described by the dependence on P, P' and D of the generalized Green's 
operator. 

We now consider the evolution of the spin amplitudes during the transport. In the simple case 
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where the time evolution of the spin does not depend on position or time, it is described by 



G S p in (t', t) = exp 



(12) 



where A spin is the evolution operator of the spin. If there is no anisotropy associated with the spin 
and if the spin is conserved during transport, then the terms (s'\ G sp i n \s) form a (25 + 1) x (25+1) 
diagonal unitary matrix. Effects depending on the spin, like circular birefringence or dichroism, 
modify the diagonal terms of (s'\ G sp in \s) such that it is not necessarily unitary anymore. Consider 
for instance a medium with absorption length kJ 1 and index n s for the spin s. In this case, we 
have 

(r', P', s'| G spin (t', t) |r, P, s) = e -^-^ s )(t' -t) ^ _ r) ^ g) (13) 

In this expression, the index n s denotes the ratio of wave celerity for a spin eigenstate s compared 
to the celerity contained in Q. Dichroism appears when the values of k s depend on s, birefringence 
when the values of n s depend on s. It is also possible to introduce spin-flips. 

To summarize our construction of the Green's operator of the homogeneous medium we have 
the following formula for Go 

(r', P', s'| G (t', t) |r, P, s) = £(r'-r, t'-t) | G spin (t', t)\s) 5 ± (DP- 1 )5 ± (p- 1 P / ) e^'"^. 

(14) 

The dependence of the Green's function as a function of the distance is, according to the ex- 
pression (fT4l) . the product of the scalar Green's function Q and the Green's function of the spin 
G sp i n . If there is no absorption, the Green's function follows the well known r _1 decrease of the 
Green's function in three dimensions. In this situation, the enhanced Green's function (TT41) verifies 
the conservation of energy by construction. 

IV. DESCRIPTION OF A SCATTERING EVENT 

In the previous section we have constructed the Green's function for a homogeneous medium. 
In this section, we include the description of the scatterers into the same formalism. The problem 



of scattering with spin particles was originally discussed by Jacob and Wick [35] when two parti- 
cles with spins are colliding. We consider more general scatterers in our formalism and describe an 
arbitrary interaction in the far-field. In this situation, we only need to know the so-called on-shell 
"T-matrix", or transition matrix, in a theory of multiple scattering. 
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We use the following definition of the on-shell T-matrix, in the case of light scattering: Con- 
sider an incoming field E in and a scatterer at position r. At a far-field position, the scattered 
field E sca depends linearly on E in 

E sca (r', = G(r', t'; r, t)TE in (r, t). (15) 

The T-matrix can be computed or measured in several ways that we will not discuss. In general, 
the T-matrix depends on the incoming and outgoing directions and on the polarization. These 
dependences will be rigorously taken into account in the present formalism. 

The position in space of the scatterers are denoted by iv A scattering event corresponds to 
a change in the direction of propagation at from P„z to P n+ iz. In the laboratory frame, the 
rotation matrix corresponding to a scattering event is R = Pn+iP" 1 . However, it is the local 
rotation which is physically relevant in a local description of scattering. The expression of this 
local rotation is R = P~ x RoP n = P~ 1 P n+ i. This is the local rotation experienced from the point 
of view of the particle during the scattering event. Using the decomposition of R into Euler angles 
(</>, 0, ip) (Eq. (OQ)), we rewrite P n+1 = P n R as 

P^Z&r 1 =P n Z(4>)Y(6). (16) 

9 is the scattering angle and <p and — tp will act as spin angles in ©. As a consequence of this 
relation, the modifications of the spin angles of P n and P n+i transform according to 

P„ -> P n Z(0) (17) 
P n+1 -> P„+iZ(-<0) (18) 

and cast the rotation into local frames adapted to the scattering event (see figure 12). 

If there is a spherical scatterer at position r, the scattering operator between incoming and 
outgoing states is 

(r', P n+1 , s'\ f (f , t) |r, P n , s) = 6®(t> - r)5(t' - t)e is ^ (Y(9), s'\f |I, s) (19) 

(Y is a rotation around the Y-axis of the reference frame) We have used the spherical symmetry 
of the scatterer to bring the incoming and outgoing frames back into the incoming one, so that the 
outgoing frame is simply Y(0). In this formulation, (Y(0), s'\T |I, s) is just the usual T-matrix 
written in the spin eigenstates basis, which we write T ss >(u, 9). We impose that there is a scatterer 
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FIG. 2: A single scattering event. The freedom of rotation of the incoming and outgoing frames is compen- 
sated by a local rotation of the spin angle ©. The azimuth angle <fi of P^Pn+i determines the scattering 
plane and the zenith angle 9 of P~ 1 P n+ i is the scattering angle. 



at position r with the density operator (r'| g |r) = <5^ 3 ^(r' — r)g(r), with 



N 



g(r) = Y,^ 3) (r-r l ): (20) 



If there is no scatterer at position r the scattered field vanishes. The scattered field is obtained 
using the the closure relation © 

(r', P', s'\G (t', t)fg\r, P, s) 

= [ dV f dP" V (r', P', s'\ G (t, t) \r", P", s") (r", P", S "| fg |r, P, s) 

JV JSO(3) s "=^s 

s 

= g(r'-r,t'-t)5 ± (DP'- 1 )g(r) ^ (s'\ G 8pin (t' - t) \s") ^ s '- s '>' j s ^% s ,,{uj). 

s"=-S 

(21) 

To obtain the last equality, we have used the fact that the spin angle of the rotation matrix P" is 
arbitrary and does not play any role in the final expression and we have therefore used the spin 
angle of P' to simplify the formula. The angles (</>, 9, ip) are the Euler angles of P _1 P'. 

The T operator in formula (T2TI) stands for a single spherical scatterer, correlations or the influ- 
ence of the other scatterers are, so far, ignored. The imaginary part of T stands for extinction. The 
extinction cross-section for spin s is then 

a s ext = --lmT ss ^,9 = 0). (22) 

UJ 
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It is constructive to compare er ext with the scattering cross section for spin s 

g 

< a « = \j2 r^9d9\T ss/ ^,9)\ 2 . (23) 
s / = _ s Jo 

The total extinction and scattering cross sections are obtained by summing over the spin eigen- 
states: a ext = ^ s °"ext and °"scatt = Z^scatt- B y conservation of energy the ratio a = cr SC!ltt /a ext 
cannot exceed 1. The ratio a is called the albedo. If a ext = a scatt , all energy captured by the 
particle is scattered. This is known as the optical theorem. Let us remark that the inequality a < 1 
valid for the total albedo is not necessarily true for the albedoes a s = er s s catt /erf xt at fixed spin. 



V. THE BERRY PHASE 



At this point of the discussion, it is interesting to point out that our description of the tra- 
jectory is able to keep track of the Berry phase of the wave. Originally, the Berry phase was 
proposed as the phase factor that can appear after a cycle during an adiabatic time evolution of 
a non-degenerated quantum state [|29|1 . Later it was discovered that the concept applies to light 
polarization as well 1 31 , 3Jj| and that the time evolution of the system does not need to be cyclic or 
unitar y 07I1 . More recently, it was shown that the Berry phase also shows up in multiply scattered 
light 04I1 . The expression for the Berry phase $ involves two factors: The spin s of the particle 
and a solid angle (which we call the geometric phase), or more generally the curvature of the 
phase space of the system enclosed by the evolution of the system during one cycle. The Berry 
phase expression is [29] 

$ = -sfi. (24) 

In polymer physics, time is replaced by a curvilinear space coordinate and the geometric phase £1 is 
called writhe |p8|l . The fluctuations of the writhe induced by thermal fluctuations of the polymer's 
shape have been studied numerically using a Monte-Carlo approach [39]. For long polymers, the 
distribution of the Berry phase was found Gaussian. There is no exact approach for the statistics 
of the Berry phase for random paths. 

For simplicity, we consider the equivalent of a nearly adiabatic time evolution and compare the 
phases of the field before and after a path for which the initial and final frames are equal. Such a 
path corresponds to a closed circuit in the phase space pill . We also assume that the scattering is 
elastic and that the spin is conserved in both scattering and propagation. Between two scattering 
events, the state of a particle is simply defined by the direction of its momentum, and the phase 
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FIG. 3: Example in the case n = 6 of the Berry phase for a trajectory with identical initial and final 
frames. (Top) Trajectory in real space. The thickness of the line changes along the trajectory to give a three 
dimensional effect: The closer it is to the observer, the thicker the line is. On each segment is represented a 
frame P, such that PjZ is the direction of the segment. (Bottom) The path is represented on the unit sphere. 
Each segment with constant direction PjZ is a point on the sphere. Two succesive points are connected by 
an arc of a circle. The area enclosed by the path on the sphere is the geometric phase Q, (in grey). The 
angles 0j (see text) have been chosen exaggerately large for the readability of both parts of the figure. In 
strong forward scattering, the trajectory on the top part is very close to a straight line and the area enclosed 
on the sphere of the bottom part is concentrated around Poz. 

space is thus the unit sphere. The adiabatic evolution would corresponds to a continuous and slow 
movement of the direction of the momentum on the sphere. In case of strong forward scattering 
each scattering event corresponds to a small change in the direction of propagation, which nearly 
satisfies the adiabatic condition for the Berry phase to occur. The geometric phase is the solid 
angle enclosed by the unit wave vector p = Pz along its trajectory on the unit sphere. 

To reveal the Berry phase, we consider a trajectory made of n scattering events, and the n+ 1 lo- 
cal frames associated to each propagation. We denote by ((pi, 9i, ipi) the Euler angles of each frame 
after the decomposition of equation (O. An example of this construction, with n = 6, is shown 
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on Figure[3] We consider the case of strong forward scattering so that the coordinates (6> i+1 , 
and (6i, (pi) of the frames P i+1 and P, are close to one another. In Figure [31 it corresponds to the 
case where the points PjZ and Pj+iz on the unit sphere are close to each other, compared to the 
radius of the sphere. We call (fa, 9i, ipi) the Euler angles of the rotations P^~ Pj + i in the same way 
as we have done in Eq. (fl~6l) . The strong scattering assumption corresponds to 9i <C 1. To ensure 
readability of the Figure [3l we have represented a path with large angles 9i that do not correspond 
to strong forward scattering. By analogy with a spinning top, for which the instantaneous rotation 
around the third axis is expressed by <p cos 9 + ip [40, chap. VI] we deduce that the rotation around 
its direction of propagation experienced by the particle is fa + fa ~ (fa+i — fa) cos 9i + ip i+ i — ipi. 
The Berry phase for a given component of the field is, in our case, identified as the phase differ- 
ence between the values of this component expressed in the two frames Po and P„. The original 
definition of the Berry phase demands that the evolution be cyclic. So we have to impose that 
the frames P and P n are equal, which corresponds to one cycle. Therefore we have <p n = fa 
mod 2ir and ip n = ip mod 2n. The extra phase for a spin state s between the frames and n is, 
according to the phase factor e 1 ( s1 l>+ s '<t>) in equation (|2Tj) . 

n-1 

$ = s y %2<f> i + $ i . (25) 

i=0 

We have used the conservation of the spin (s' = s). We rewrite the sum using the quantities 
(1 — cos 9i) and we get, up to a multiple of 2n: 

n—l 

$ = -sj^(l - cos9i)(fa +1 - fa) = -sQ mod 2vr. (26) 

8=0 

We recognize that the sum in (l26l) is a spherical area. It is equal to the geometric phase up to 
a multiple of 47r because of the indetermination of the solid angle [39]. We have retrieved the 
expression of the Berry phase (|24|) . The equality between the sum in expression (|26l) and the geo- 
metric phase holds modulo An, rather than 2ir, for topological reasons 112 811 . The indetermination 
modulo 47r plays here no role because s is an integer or a half-integer, so the phase is undetermined 
modulo 2n. 

Remarkably, the Berry phase emerges in our local frame model without any special effort. The 
geometrical nature of the Berry phase appears here clearly. As it was already demonstrated, the 
existence of the Berry phase is not restricted to cyclic systems and can be extended to non close 
path and to non-unitary evolution (in our theory, if there is absorption, for instance) 113 VTI . The 
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42J. Our 



statistics of the Berry phase for random walks is a technical and difficult problem H41L 
approach will make it possible to obtain exact results concerning the Berry phase statistics. 

A possible experiment to observe the Berry phase of light consists in setting a directional lin- 
eary polarized source at the edge of a slab of a medium filled with scatterers and observe the 
outgoing light through a collimator and a linear polarizer. To ensure strong forward scattering the 
scatterers must be large such that the anisotropy g, defined as the average cosine of the scattering 
angle, verifies 1 — g <C 1. One can for instance place the polarizer in front of the source, its plane 
orthogonal to the direct beam and its axis making an angle 9 with the polarization of the source. 
The intensity as a function of 9 will display deviations from the expected cos 2 9 behavior if there 
were no scatterers. If the slab depth is large, the Berry phase will spread out on a width larger 
than 2n and the intensity will be constant as a function of 9. The experiment has therefore to be 
performed with slab depths of a few transport mean free paths £* = £/(l — g) (£ is defined in the 
next section). Different concentrations will lead to different deviations. 



VI. GENERALIZED GREEN-DYSON EQUATION 

The complexity in multiple scattering stems from the summation of the contributions to the 
field of all the paths from the source to the observer, and potentially from their interferences. In 
a first approximation one often neglects these interferences and assumes that the scattering events 
are independent from each other. In a disordered system the averaging over disorder can thus be 
performed without the sophisticated diagram techniques needed to preserve interference effects. 
The elements of the theory that we have introduced in sections [III] and [IV] are used in this section 
to write a generalized Green-Dyson equation for the multiple scattering of polarized waves under 
the assumption of independent scattering events. This equation applies to directional Green's 
operators: It relates the free Green's operator Go we have introduced in section [HI] to the effective 
Green's operator G of the same medium filled with scatterers. In the second part of this section, 
we introduce the Fourier transform of the generalized Green-Dyson equation. 

Consider a medium filled with scatterers at fixed positions The presence of the scatterers 
is equivalent to the introduction of the perturbation in the transport equation. In the formalism of 
Green's functions it corresponds to introducing a transition matrix, or T-matrix. We now consider 
the perturbed Green's operator G developped as a Born expansion in terms of Go: 

G = G + G t^Go + G T£G T£Go + ■ • ■ (27) 
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The first term of the right-hand side of (1271) stands for the unscattered field, the second term for 
the single scattering and so on. 

For a known distribution of scatterers, the operator G is the Green's function of the multiply 
scattering system for polarized waves. If the functions Q, T ss (cu, 6) (defined in section ITU and HVl 
respectively) and the operator A sp i n are known, equation (ITTT) for G solves the problem. But usu- 
ally we are interested in disordered systems, for which the exact position of scatterers is unknown. 
It will thus be impossible to make a prediction depending on the particular realization of disorder, 
but it is possible to average over all such realizations (ensemble average) and to write a general- 
ized Green-Dyson equation for the average operators. We formulate the results using the simplest 
possible model of disorder. To perform the average over disorder, we neglect the volume of the 
scatterers and consider that each position r-j is uniformly distributed in the volume V. The contri- 
bution of disorder in (1271) comes exclusively from the averaging of the operator Tg. We introduce 
the self-energy S(o?) 

EH = T(u)q(t). (28) 

The line denotes the average over all positions for all the scatterers, which we will explicitely 
define now for the calculation of the average density. The average of an operator O is simply 



defined for arbitrary (1| and |2) by (1| O |2) = (1| 6 |2). If the scatterers all have the same 
operator T(ui) the self-energy becomes 

±{u) = T(^W = t(u,) J v ^-J v yf^^-r) = TH f; I = f ee pTM, 

(29) 

where p is the mean density. Our approach of disorder is nonetheless naive, we have neglected 
the possibility that a trajectory visits several times the same scatterer and thereby omitted some 
possible correlations in the averaging procedure. The topic of correlations created by several visits 
to the same scatterer is addressed abundantly in the literature. The most remarkable result is that 
for a sufficiently small density of scatterers, the expression of the Green-Dyson equation remains 
the same, but with a modified self-energy expressed as an expansion in powers of p [7]. The 
first order of the expansion is exactly the expression we have found in formula (1281 ), which is 
an encouraging fact in favor of an equivalent expansion in our framework. Such an equivalent 
expansion would extend the domain of validity of our results to higher densities. We will address 
this question in our future work. 

Let us call G the average Green's operator resulting from the averaging over disorder of the 
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operator G defined in (|T7l) . In the case where all scattering events are independent, the average of 
the products involving several times £ are the products of the averages, this is sometimes called the 
Twersky approximation The Born expansion (|27T) becomes, in this approximation the series 

G = G + GqSGq + GoSGqSGo + • • ■ (30) 

We have so far expressed the functions only in the direct space of position and frames. Physical 
problems however are easier to formulate in the reciprocal space, because convolutions become 
products. We denote by |q) the reciprocal space representation of |r) and we have the well known 
contraction 

(r|q) = e iqr (31) 
which gives, together with the closure relation J d 3 r |r) (r| = 1, the fundamental relations for the 



Fourier transform 04711 . Similarly, we use the representation in reciprocal time that has the same 
properties as the reciprocal space. 

As the system is invariant along time translations, the operators T, G , G, noted as 6, depend 
on one frequency u. Similarly, in the average medium obtained after integration over disorder, the 
translational invariance is restored the operators only depend on one vector q. We introduce the 
notation 0(u>, q) for operators acting on the Hilbert space made of the elements |P, s), so that 
(q| 0(w) |q) = 0(uj, q) ^ 3 )(q' — q). The representation in reciprocal space for rotation and spin 
has different properties and is presented in the next section using an expression equivalent to (|3TT ). 

The Born expansion written using the operators £(a;), Go(q, uj) and G(q, u>) reduces to 

G(q, lo) = G (q, u) + G (q, w)E(w)G (q, u) + G (q, w)£(a>)G (q, cj)E(o;)Go(q, to) H 

(32) 

The operator £(cg>)G (^, q) has a norm ~ pcr ext g _1 and is therefore smaller than 1 for large q or 
sufficiently small density p. Consequently, the series 5^ n (E(a;)G (c<j, q)) n converges and we can 
rewrite the expansion (1321) in a self-consistent way 

G(q, cj) = G (q, uj) + G (q, ou)±(u)G(q, u). (33) 

Equation (|33l is a generalized Green-Dyson equation for the directional Green's operators. The 
convergence of the series occurs when qi > 1 where i = \j pa ext is the extinction mean free path. 
The theory described by the generalized Dyson equation does not apply to situations where q£ < 1 
because scatterers are not in the far-field of each other anymore. In the next section, we introduce 
the rotational harmonics transform on the variable P to solve equation (1331) . 
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In the expression of the free Green's operator Go(u, q), Eq. CHI), the factor 5j_(DP _1 ) mod- 
ifies the usual expression of the Green's function ~ l/(^ — q 2 ± rr/), where the quantity 77 is a 
infinitesimaly small positive quantity. Using the identity 

/ iwr/c 1 1 

d 3 r^ e-^MDP- 1 ) = rr 
4vrr ±K 1 47r(^_ q .pz) 2 

we get the Fourier transform 

<P', s'\ G«> A (u, q) |P, s) = ^ (s' G spin M s) (P'|P> . (34) 

[- — q ■ PzJ ± 17] N ' 

The R, A superscript stands for the retarded and advanced Green's operator. As we have already 
mentioned, the scalar theory must be retrieved if we sum our expressions over P and s. Here we 
can check that equation (l34l) integrated over P and summed up on s yields formally 



£ / d P 1 „.v 2 ^(q. «0 = %(ci, w). (35) 

s t^ s Jso(3) (|_q.Pz) ^r-q 2 

The regularization by the infinitesimal imaginary term ±w? can be performed on the r.h.s term of 
Eq. (1351) so that we retrieve the usual Green's function for the field component with spin s = 0. The 
superposition of all partial directional fields cancels out because of the phases of all the spins with 
s ^ 0. The final result depends on q = |q| : the directivity of the directional Green's operator (l34l) 
is lost. We call the expression (- — q ■ Pz) -2 the directional factor of the Green's operator. 

We can perform the same integration as in formula (|3~51) on each operator G (u, q) 
in formula (|33l) . We retrieve the scalar Green-Dyson equation G(u, q) = Q{oj, q) + 
G(co, q)T,(uj)G(u, q) with 

£M = ip rsny6d9%o(u, 9). 
* Jo 

The Green's operator G does not necessarily act separately on |r, t) and |P, s). In optics for 
instance the Faraday effect in a medium with Verdet constant V creates a phase shift sVH ■ p. In 
this case, provided that the spin is conserved during propagation, equation (1341) is modified and 
the Green's operator Go becomes 

(P', s'\ G*' A (q, u) |P, s) = - 2 (s'\ G spin (w) \s) (P'|P) . (36) 

+ (sVB - q) • Pz) ± irj 

The separation of the article into two parts to explain the harmonic transforms is intended 
only to keep the explanations simple. The Fourier transform presented in this section and the 
rotational harmonics transform that we introduce in the next section could have been presented 
simultaneously. This is actually mandatory in the case of a spin-orbit coupling involving |P, s) 
and q like the Faraday effect, as illustrated by formula (1361) . 
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VII. ROTATIONAL HARMONICS REPRESENTATION 

The linear representation of compact Lie groups provides an harmonical analysis on the group 
SO (3) which is the equivalent of Fourier series for periodic functions. These harmonical analysis 
is different from the Fourier transform because SO (3) is a compact, non-commutative Lie group. 
For our purpose there is no need to know this mathematical theory in detail. Some formulae used 
in our theory are presented in the appendix lAl 

The reciprocal representations are labelled by a main index I > which can be an integer or a 
half integer and two sub-indices m and n taking the values —I, — I + 1, . . . , I. We denote them by 

I mn ) and we have 



<R| In > = i" V2TTT ef* < n (cos 9) (37) 
with 0, 9 and ij) the Euler angles of R. The closure relation, 

CO I I 

EE El-)(^l = l, (38) 
l=A m=—l n=—l 

is useful to obtain most of the result provided in the present work. We call the function (R| J n ) 
the rotational harmonics, (the misleading terminology "generalized spherical harmonics" is often 
used). These functions are normalized Wigner -D-functions ll 1 8l 14311 . The indices of the sums on I, 
m and n in (l38l) have steps equal to 1. A, the lowest order of /, depends on the spin S of the wave 
and takes the value or 1/2 if the spin S is an integer or a half-integer respectively. It follows 
that the indices /, m and n are either all integers or all half-integers depending on S. For a given I 
(the order of the harmonics) there are (21 + l) 2 orthogonal functions, which seems a lot. We will 
see that for / > S only (2S + 1) (2Z + 1) harmonics components have to be taken into account. 
The functions d! m n are defined from ZYZ Euler angles and are real-valued, their full expression 
is given in Equation (IA1|) . We have added the unconventional factor i n_m to obey the identity 
(R| ^ m )* = (R -1 1 J m ) (that is we have suppressed the factor (— l) n_m in this latter relation). 

The rotational dependence of a state can be expressed as a function of the frame P with the 
bracket notation (P | or as a function of rotational harmonics indices ( ^ n | v|/). Our computation 
uses this latter representation because it is related to the spin in a fundamental way: Let us use the 
closure relation of the rotation representation J so ^ dR |R) (R| = 1 to get the expression 



i n ,s|vl/>= / dR(^,s]R, s)(R, s|V) 

JSO(3) 



(39) 
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According to the relations © and (1371) the integration over the spin angle just leaves a factor 8 ns , 
stating that the representation for a spin eigenvalue s involves only the kets \ r j ls ). In the theory 
of multiple scattering, S is fixed and there is only one species of particle involved. It means that 
in the Hilbert space of physical states the equality | r l nn ) ® \s) = | J^) <S> \s) 5 ns is valid and the 
n-index can be left out without any restriction on the physical content of the theory. 

We use the rotational harmonics to expressed the products of the operators used in the general- 
ized Dyson equation (1331) : 



K. \G R ' A {q, 
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Would we have used the representation in SO (3), we would have obtained integrals instead of the 
discrete sums. The sums over L and L' are infinite. However, like in Fourier series, it is possible 
to truncate the sum because the brackets amplitudes decrease for large L like a L (for a certain a, 
< a < 1) such that high order coefficients can be neglected. The sums on indices L, M and 
a can be seen as sums over a single index. Equation (|40l) is therefore a discrete linear equation 
and is, up to a one-to-one transformation of (L, M, a) into a single index, a matrix equation. 
If the coefficients of Go(q, lo) and E(a>) are known, one can compute G(q, w) with standard 
linear algebra and get a general solution for the transport of polarization in multiple scattering by 
inverting the matrix of 1 — G E 

The directional Green's operator G describes the transport in a effective, homogeneous 
medium. Paths statistics and Berry phases are included in this description, which represent an 
improvement over the usual effective Green's functions used in multiple scattering theories. In the 
presence of several different kinds of scatterers, with densities and scattering matrices Tj, we get 
an effective Green's operator by replacing E by piTf. The scattering properties of the different 
scatterers are averaged and their respective weights are proportional to their respective densities. If 
the scatterers are not spherical and have random independent orientations, the effective T-matrix 
is the average over all orientations of the orientation-dependent T-matrices. In particular, if the 
orientation probability is uniform, the average T-matrix is equivalent to a spherical scatterer's 
T-matrix and the spherical symmetry can be used as is described in the next section. 

To compute the rotational harmonics expansion of G up to order L, one needs to use approxi- 
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mately (2S + 1)L 2 triplets (I, m, n) for the incoming frame, and as much for the outgoing frame. 
Therefore required number of rotational harmonics coefficients for the computation of the Green's 
operator scales like (2S + 1) 2 L 4 . The complexity of the inversion of the linear system (|40l) scales 
approximately like L 5 . 

VIII. THE ROLE OF ROTATIONAL INVARIANCE 

We have considered in the previous sections the general case of an homogeneous medium 
invariant under translations and we have decomposed the vectors of the three-dimensional space 
into a radius part and an angular part to construct the general formalism for a directional Green's 
operator. A special and important case shows up when the medium and the scatterers are also both 
invariant under rotation. In the SO (3) representation, we have in this case for all frames P and P' 
the equality already discussed in section [IV] 

(P', s'\ t(u) |P, s) = (P^P', s'l |I, s) . (41) 

Using the closure relation (1381) we get 

oo I 

<p', s'l eh ip, s) = e E ( p ' lp '^'\±)(±\ t ^)\ 1 ^)^ 

Z=A m,n=—l 

We expand the (P _1 P' | ^ n ) thanks to formula (|A11|) and sum over n. After using formula (IA12|) . 
we obtain 

oo I 1 

(P', 5'| t(u) |P, s) = E 7|f=y< p, | P i'><w|SM|/, S )<P|A>. 

i=A m,p=—l 

By identification of the terms of the right-hand sum with the rotational harmonics expansion in P 
and P' of the left-hand term (IA7I) we get 

/ ;' I v^/ \ i ; \ IV Omm 1 / ; | A/ \ t T v ^ll'^mm' / / \ /a^\ 

< «' 1 EM 1 ~> = vim ( 1 E( "» s> s VITT 2 A M ' <42) 

If the medium and the scatterers are both invariant under rotations, the free Green's operator G 
and the Green-Dyson operator G are both independent from the reference frame. These operators 
depend on q and on the incoming and outgoing directions of propagation Pz and P'z, a rotation 
of the reference frame acts simultaneously on these three vectors. Using the decomposition 

q = <?Qz, (43) 
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where Q is a rotation matrix and q = |q|, we express the invariance of the Green's operators as 



P', 8 



<R,Al 



Q- A P,s>, (44) 



where G stands either for Go or for G. In Equation (l44l) . the operator G depends only on q and u. 
Using formula (1341) . one can even simplify further the Green's operators by defining 

2 

G ' [q, w) = — go ( — ) (45) 



G R %, U ) = ^ f^) . (46) 



-.2 

As a consequence of formulae (fl4l) we find a symmetry property for g 



V I I I 

m's' I 60 mn / ^ "mm' 



, } OC 5 r , 



The same symmetry has been observed for the T-matrix in Equation (1421) and we conclude that it 
will propagate to g. A straightforward calculation shows that Equation (|33l) takes the form 

The non-vanising rotational harmonics coefficients are related by 

/ i' I pR,A v R > A (Sl\\ i \ + 

\W|o ^ ^ I ms I — \ ms' I 60 V Co" / 

L=A V (T= — S cr' = — 5 

Thanks to invariance under rotations, we have obtained a simplified formula for the generalized 
Green-Dyson equation. The matrices E ) are already known for several spherical scatterers, like 
Mie scatterers. Formula (l48l) contains only three sums, two of which, over the indices a and a', 
are finite sums and correspond to a matrix product. The expansion of g into rotational harmonics 
coefficients is therefore a single sum over a generalized index instead of a double sum as in (|40l ). 
The computation of g requires to know the rotational harmonics coefficients of the S-matrix and 
of g . The coefficients of g can be computed starting from expression (|34|) : 



ms' 



2 V 2 y_ x (1 -/ix) 2 ±ir/' 

For arguments x > 1 there is a pole at 1/x in the integral (|49| ) which has the same physical origin 
as the pole of the Green's function in a scalar representation. 
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IX. CONCLUSION 

We have presented a model for the transport of a polarized wave in a complex medium con- 
taining anisotropics. In our approach, the wave is represented in the particle picture of quantum 
mechanics; the transport of the wave is described as a superposition of trajectories with attached 
probabilities. The polarization of the wave follows from the spin of the particles. We express the 
amplitudes of the spin eigenstates in local frames. A local frame is a rotation matrix which third 
axis is colinear with the momentum of the particle. Local frames form a convenient setting for 
polarized waves because the phases of the spin eigenstates of a particle are related to each other 
by frame dependent factors. Based on the local frame representation, we have established a theory 
for the multiple scattering of particles with spin and we have obtained a generalized Green-Dyson 
equation which takes into account the three dimensional nature of the transport. One can real- 
ize the three-dimensional nature of the formalism by noticing that the Berry phase, by essence a 
three-dimensional concept, naturally emerges from our equations. 

We have suggested to solve the generalized Green-Dyson equation using rotational harmonics, 
which transforms the convolutions on SO (3) into matrix products. The solutions of the generalized 
Green-Dyson equation are obtained by linear algebra operations. The local frame description is 
able to take into account several circular anisotropics, like birefringence and dischroism, Faraday 
effect and anisotropic scattering. In the special and important case where the scatterers and the 
medium are both invariant under rotations, the generalized Green-Dyson equation takes a simpler 
form and should be solvable with reasonably light numerical power. Linear birefringence and 
dichroism can also be described in the theory but require a more complex mathematical treatment 
which was not presented here for the sake of simplicity. Beyond the theory of multiple scattering, 
we believe that the generalized Green-Dyson operator is a useful object for the statistical study 
of random three-dimensional paths, thanks to its ability to perform path integrals of direction 
dependent functionals. 

The computation of the solution of the generalized Green-Dyson equation using rotational har- 
monics is probably the most efficient, because it reproduces the properties of the Fourier transform 
used to solve the usual Green-Dyson equation. On one hand, as the rotational harmonics form a 
discrete basis of the algebra on SO (3), one only needs to compute for each value of cq/uo a discrete 
set of coefficients. On the other hand, the maximum order L one has to use to get a physically 
relevant result grows with the strongest anisotropy of the system. Despite the simplicity of the 
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expression of the solutions, it is not conceivable to draw the calculations analytically in general. 
The numerical difficulties lie in two points: The inversion of the linear system ([40]) or (1481) and the 
expression of the operator g (l49l) . 

We have neglected the interactions between scatterers, which can create correlations between 
their positions. Interactions can be handled by the introduction of the pair correlation function (or a 
structure factor) which would enter as a factor in the second order term of the Born expansion[44]. 
In the terms of higher order n of the Born expansion, one should manipulate n-points correlation 
functions, which introduce a supplementary difficulty. Much less is known concerning the inter- 
acting scatterers than concerning the self-energy. We also considered that the medium is globally 
invariant under translation and that the density of scatterers is homogeneous, which excludes small 
size systems and leads out the effects of boundaries. Our theory concerns the bulk of a system, but 
experiments are usually carried out from the edges. It is therefore an important issue to investigate 
the role of the boundaries. 

As such, our work should rather be considered as a first step towards a full theory of polarization 
transport in presence of anisotropics. It may nonetheless already lead to experimental applications 
and investigations in complex media. On a more general point of view, our work suggests that 
polarized waves in complex media contain more information concerning anisotropics than scalar 
waves and that polarization is therefore a relevant supplementary observable for the investigation 
of anisotropic systems. Our future work will be dedicated to improve the domain of application of 
the presented framework to denser, interacting or quantum systems. First of all, the expression of 
the self-energy expansion in powers of the density should be extended up to higher orders and the 
formalism for intensity diagrams should be developed. 
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APPENDIX A: THE ROTATIONAL HARMONICS 



The rotational harmonics analysis is similar to the Fourier analysis for periodic functions. It 



is the direct application of the theorem of Peter and Weyl |45|] to the representations of the group 
SO (3). Peter- Weyl theorem states that any complex function / defined on SO (3) may be devel- 
opped into rotational harmonics series. The representations of SO (3) were introduced in physics 



by Wigner fl 1 8M but remain confidential in physical theories mostly because in spinless theories the 



rotational harmonics simplify to spherical harmonics. A important application in physics was to 
establish Wigner- Eckart's theorem Jisil^] . (One can notice that the equation (|42|) is a particular 
case Wigner-Eckart's theorem). It was shown in section [VTII that the set of rotational harmonics 
useful in the problem of scattering of spin particles depends on the spin S. More precisely, it 
depends whether S is an integer or a half odd integer. The rotational harmonics with a half odd 
integer main index stand for fermions while the ones with an integer main index stand for bosons. 
The rotational harmonics are given by formula (|37l) and 



min{£+m, I— n} , g\2l+m—n—2p i . g\n— m+2p 

d l mn (cos9) = - m)\ (I + m)\ (I - n)\ (I + n)\ V (-If V cos 2^ 

p=max{0, m—n} 

(Al) 

Remark that the arguments of factorials and the powers are integers, in both cases where I, m and 
n are integers or half-integers. The computations of the set of d l mn functions is simplified by the 
symmetries 

< n (cos0) = (-l) l+m d l mj _ n (-co S 9) = (-l) i+m <_ n (cos(7r - 6)), (A2) 
= (-l) m -V mi _ n (cos^), (A3) 

= (-ir-X m (cos0). (A4) 

One can also define the rotational harmonics from the spherical harmonics by mean of an operator 
introduced in references [19]. This result is sketched in the next appendix. 
In our notations the harmonics expansion is written 

i 

/(R) = (R|/> = E E ( R \±)(±\f)> (A5) 

/>A m,n=—l 

where the harmonics coefficients are 



' i 

i mn 



\f) = [ dR (^JR) (R|/> = f dR (R\±y /(R). (A6) 

JSO(3) JSO(3) 



27 



Similarly, the operators can be expressed in the basis of the rotational harmonics. 

0(R',R) = <R'|6|R> = E E (Jn'\0\^ n ){R'\J nl )(R\^ n r, (A7) 



l,l'>A m,n=—l m'n'=—V 



where < J n , \ 6\±)= [ dR f dR 0(R, R) <R| J n , }* (R|j n >. (A8) 

-/,SO(3) J50(3) 

The product of operators C = AB is expressed in direct space as a convolution 

(R| C |R) = / (R| A |X) (X| B |R) dX (A9) 

J SO (3) 

and in harmonics space as a series 

oo L 

( m'n' I ^ I mn ) = E/ E/ ( mV I ^ I JWf ) ( MN I ^ I mn ) • (A10) 

L=A M,N=-L 

The sum over M and N is the usual matrix product. In our calculations we also use the relations 

J2 (Ri I 4 > <R 2 \ P l n ) = V21 + 1 (RxR 2 1 ± ) , (All) 

p 

<R- 1 |^ n ) = (R|^)*. (A12) 
It is also interesting to note that we have the Parseval-Plancherel formulae: 

„ oo I 

(f\f}= / |/(R)| 2 dR = ^ Yl K"inl/)| 2 > (A13) 

Js °(V l=A m,n=-l 

/dR/dR |(R|6|R)| 2 = £ E \(Jn>\0\±) 

l,l'=Am,n=—lm',n'=—l' 



(A14) 



APPENDIX B: ROTATIONAL HARMONICS UNDER OTHER NAMES 

In this appendix we compile the relations between the rotational harmonics and other spe- 
cial functions used in literature under several names and with a large variety of normalizations. 
This appendix is intended to allow those familiar with one the numerous forms of the rotational 
harmonics to comprehend our presentations from a wider point of view. Conventionnally, the 
argument cos 9 may be written fx. 

Wigner introduced the rotational harmonics arranged in matrices called Wigner's D-matrices 



of order I [|18ll . The elements of D(l) are given by 

D l m , m {<j>, 9, V) = e-^^JcoB^e-*™*. (Bl) 
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FIG. 4: Examples of functions d l mn (cos 9) for small values of I. Left: graphs of functions with spin index 
for I = 1 and I = 2. One recognizes the unnormalized spherical harmonics. The only function with I = 0, 
dg = 1 is not shown. Middle: graphs of functions with spin index 1 for I = 1 and I = 2. Right: graphs of 
functions with spin index 1/2 for I = 1/2 and I = 3/2. 

1. Functions tt and r used in the Mie expansion 

The well known Mie theory for spherical dielectric scatterers provides the general solution of 
the scattering of a plane wave by a dielectric sphere of arbitrary size and refraction indexyj, |2fl. 
Forward scattering appears for large spheres of radius a ^> 2n/X, this is sometimes called the Mie 
effect. The scatterers are spherical, so they are invariant under rotations and the T-matrix is thus of 
the form of Equation (|42l) . Written in the spin eigenstates basis (also called circular polarization 
basis), the T-matrix coefficients take the form 

T ss ,(u, 9) = i s (Sx(9) + ss'S 2 (9)) 8, s 1 = ±1 (B2) 

u 

with the coefficients Si(9) and S 2 (9) as defined by van de Hulst |4j]. 

Mie theory provides the expansion of the coefficients Si and S 2 into series of terms involving 
special functions referred as n n and r n in Q4£| and many other textbooks. The coefficients of the 
expansion are called a n for S± and b n for 5*2. In the circular polarization basis, we obtain the main 
coefficients as combination of n n and r n : 

S 1 (9) + S 2 (9) = E^=i(2^ + l) ^'" ( ^ + + T) (CMg) » (B3) 

S 1 (9)-S 2 (9) = J2Zl + 1) °" 2 ^ " n(C ° S n(n+l) COS ^ ■ 

The functions 7T/(//) and Ti((jl) are defined by 

7T = 0, 7Ti = 1 and 7T; +1 (/l) = 21 ^ 1 /i7T;(/i) - ~-~~7T;_i (/^) ] (B4) 

n(x) = Ihtti(ii) - (I + i)n-i(jJ>)- (B5) 
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As it was shown by Domke Q20Q, we have the relations: 

7r,0i) = 1(1 + 1) (d^M + 4,-iC")) , (B6) 
r,(/i) = 1(1 + 1) (^(/i) - • (B7) 

(In [|50, the functions n and r are defined with the same equations, without the factor 1(1 + 1).) 

It appears that the Mie expansion is nothing but the rotational harmonics expansion for spins 
s — ±1 of a spherical scatterer's T-matrix. In the low density approximation used in the article, 
the self-energy matrix defined in (l42l) is then 



£^=-£-^=2^(2/ + !)^^. 

U) 2 

One should not be surprised because Mie's result comes from the resolution of Maxwell's equa- 
tions, which are vectorial, thus by nature concerning fields of spin 5=1. 



2. Vector spherical harmonics 



In electromagnetism, the introduction of spherical harmonics is motivated by the possibility of 
multipole expansion to solve several problems Jj, ofl . The multipolar expansion of a vector field 
can be expressed in terms of "vector spherical harmonics". There does not seem to be a standard 
notation for these vectors. In Ref. yQ, one finds the definition of X/ m (0, (p) while in Ref [5] we 



have 



r(e) 
Im 



r) 



1 



--rVY lm , 



y, ( :V) 

Yg»(r) 



The spherical harmonics Yi m being given by 

Y lm (6, <P) = 



Yim r - 



f X Y 



(e) 
Im ' 



21 + 1 (l-m)\ 
4tt V (l + m)\ 



<4,o(cos0) e 1 ' 



we have in spherical coordinates (r, 6, <j>) 



YS(r) 



(-l) m i l 21 + 1 (l-m)\e 



im<f> 



y/l(l + 1) V 4tt y(/ + m)!sin# 



1=1*1 (Al 



\ 



\ 



(d l m ,M + d m,-M) 



(B8) 

(B9) 
(BIO) 



(Bll) 



(B12) 
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Only the vector spherical harmonics with m = ±1 are used to solve electromagnetic problems. 
We find again that it is more natural to look for the solutions using the basis adapted to spin 5 = 1 
in electromagnetism. 



3. Generalized spherical harmonics 

In the work of Kuscer and Ribaric, "generalized spherical harmonics" have been introduced 
for the Stokes parameters describing polarized light Stoked parameters are intensities and 
therefore correspond to the amplitude of the wave field squared. This is why the rotational har- 



21 



2311 . In other words, the 



monies appear in this work and followers with a power of two [20, 
expansion involves the functions d l m0 and d l m±2 . The spherical functions are, in these papers, 
noted P l mn (cos9) and are simply 

P l m M=i n - m d l mn {ri. (B13) 
Our rotational harmonics (R| J n ) only differ from these spherical functions by a normalization 



factor s/2l + 1 which we introduced to simplify the expression of the convolution formulae. 



4. Spin- weigh ted spherical harmonics S Y{' 



Newman and Penrose introduced functions called the "spin-weighted spherical harmonics" to 
study gravitational radiations (which have a spin 5 = 2) II 1911 . The spin weighted functions are 
obtained after derivation of the spherical harmonics by the operator 5 (read "thop") defined, for a 
spin- weight s function f s , by 

d 



Qf s 



sin ( 



d_ 1 
dO sin 6 > 



sin i 



7- 



(B14) 



The usual spherical harmonics are of spin- weight and the spin-weighted spherical harmonics of 
spin s are defined by s Yi m (9, 4>) = 
for s < 0. We have the relation 



d s Y lm for s > and a Y lm (6, <P) = \ 



(l+s)\\ V > Ilm 



-1 



21 + 1 Al I 

d mi _ s (cosi 



irruj) 



(B15) 



[1] M. Born and E. Wolf, Principles of optics (Cambridge University Press, 1999), seventh ed. 



31 



[2] G. Mie, Ann. Phys 331, 377 (1908). 

[3] J. D. Jackson, Classical electrodynamics (John Wiley and sons, 1972), 2nd ed. 

[4] H. C. Van de Hulst, Light scattering by small particles (John, Wiley and sons, New York, 1957). 

[5] R. G. Newton, Scattering theory of waves and particles (Springer, New-York, 1982). 

[6] M. I. Mischchenko, L. D. Travis, and A. A. Lacis, Multiple scattering of light by particles (Cambridge 
University Press, 2006). 

[7] M. C. W. van Rossum and T. M. Nieuwenhuizen, Rev. Mod. Phys. 71, 313 (1999). 

[8] E. Akkermans and G. Montambaux, Mesoscopic Physics of Electrons and Photons (Cambridge Uni- 
versity Press, 2006). 

[9] F. A. Pinheiro and B. A. van Tiggelen, J. Opt. Soc. Am. A 20, 99 (2003). 
[10] G. Labeyrie, F. de Tomasi, J.-C. Bernard, C. A. Muller, C. Miniatura, and R. Kaiser, Phys. Rev. Lett. 
83, 5266 (1999). 

[11] T. Jonckheere, C. A. Muller, R. Kaiser, C. Miniatura, and D. Delande, Phys. Rev. Lett. 85, 4269 
(2000). 

[12] D. V. Kupriyanov, I. M. Sokolov, C. I. Sukenik, and M. D. Havey, Laser Phys. Lett. 3, 223 (2006). 
[13] I. M. Sokolov, M. D. Kupriyanova, D. V. Kupriyanov, and M. D. Havey, Phys. Rev. A 79, 053405 
(2009). 

[14] G. Labeyrie, D. Delande, C. A. Muller, C. Miniatura, and R. Kaiser, Optics. Comm. 243, 157 (2004). 

[15] L. Margerin, M. Campillo, and B. van Tiggelen, J. Geophys. Res. 105, 7873 (2000). 

[16] H. Sato, Geophys. J. Int. 121, 523 (1995). 

[17] I. Kuscer and M. Ribaric, J. Mod. Opt. 6, 42 (1959). 

[18] E. P. Wigner, Group Theory and its Application to the Quantum Mechanics of Atomic Spectra (Aca- 
demic Press, New York, 1959). 
[19] E. T. Newman and R. Penrose, J. Math. Phys. 7, 863 (1966). 
[20] H. Domke, Z. Meteo. 25, 357 (1975). 
[21] C. E. Siewert, Astron. Astrophys. 109, 195 (1982). 

[22] D. Lacoste, B. A. van Tiggelen, G. L. J. A. Rikken, and A. Sparenberg, J. Opt. Soc. Am. A 15, 1636 
(1998). 

[23] J. W. Hovenier and C. V. M. van der Mee, Astron. Astrophys. 128, 1 (1983). 

[24] R. D. M. Garcia and C. E. Siewert, J. Quant. Spectrosc. Radiat. Transfer 36, 401 (1986). 

[25] S. R. Pal and A. I. Carswell, Appl. Opt. 24, 3464 (1985). 



32 



[26] A. H. Hielscher, A. A. Eick, J. R. Mourant, D. Shen, J. P. Freyer, and I. J. Bigio, Opt. Expr. 1, 441 
(1997). 

[27] B. D. Cameron, M. J. Rakovic, M. Mehrubeoglu, G. W. Kattawar, S. Rastegar, L. V. Wang, and G. L. 

Cote, Opt. Lett. 23, 485 (1998). 
[28] V. Rossetto and A. C. Maggs, Eur. Phys. J. B 118, 323 (2002). 
[29] M. V. Berry, Proc. R. Soc. Lond. A 392, 45 (1984). 

[30] A. Shapere and F. Wilczek, eds., Geometric phases in physics (World Scientific, Singapour, 1989). 

[31] A. Tomita and R. Y. Chiao, Phys. Rev. Lett. 57, 937 (1986). 

[32] M. V. Berry, J. Mod. Opt. 34, 1401 (1987). 

[33] R. Bhandari, Physics Reports 281, 1 (1997). 

[34] A. C. Maggs and V. Rossetto, Phys. Rev. Lett. 87, 253901 (2001). 

[35] M. Jacob and G. C. Wick, Annals of Physics 7, 404 (1959). 

[36] P. G. Kwiat and R. Y. Chiao, Phys. Rev. Lett. 66, 588 (1991). 

[37] J. Samuel and R. Bhandari, Phys. Rev. Lett. 60, 2339 (1988). 

[38] F. B. Fuller, Proc. Natl. Acad. Sci. USA 68, 815 (1971). 

[39] V. Rossetto and A. C. Maggs, Phys. Rev. Lett. 88, 089801 (2002). 

[40] L. D. Landau and E. M. Lifshitz, Mechanics (Butterworth Heinemann, Oxford, 1993), 3rd ed. 

[41] M. M. G. Krishna, J. Samuel, and S. Sinha, J. Phys. A: Math. Gen. 33, 5965 (2000). 

[42] V. Rossetto, Europhys. Lett. 69, 142 (2005), arXiv:cond-mat/0409268. 

[43] I. M. Gel'fand and Z. Y. Sapiro, Am. Math. Soc. Translations 2, 207 (1956). 

[44] S. Fraden and G. Maret, Phys. Rev. Lett. 65, 512 (1990). 

[45] F. Peter and H. Weyl, Math. Annalen 97, 737 (1927). 

[46] C. Eckart, Rev. Mod. Phys. 2, 305 (1930). 

[47] There are several presentations of the Fourier transform, with different normalizations. Here, the nor- 
malisation consistent with our formulas is f d 3 q |q) (q| = (27r) 3 l. 



